Depletion-induced percolation in networks of nanorods 
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Above a certain density threshold, suspensions of rod-like colloidal particles form system-spanning 
networks. Using Monte Carlo simulations, we investigate how the depletion forces caused by spheri- 
cal particles affect these networks in isotropic suspensions of rods. Although the depletion forces are 
strongly anisotropic and favor alignment of the rods, the percolation threshold of the rods decreases 
significantly. The relative size of the effect increases with the aspect ratio of the rods. The structural 
changes induced in the suspension by the depletant are characterized in detail and the system is 
compared to an ideal fluid of freely interpenetrable rods. 
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Networks of electrically conducting carbon nanotubes 
(CNTs) have important applications in the develop- 
ment of light-weight conducting composites and films 
Q, HIE HQ- At equilibrium, colloidal suspensions of 
CNTs contain clusters with a distribution of sizes. The 
contacts within each cluster extend conductivity beyond 
the length of individual CNTs, and contiguous conduct- 
ing paths first appear on a macroscopic scale when the 
CNT density crosses the percolation threshold, at which 
the average cluster size diverges. The conductivity of the 
sample rises sharply by several orders of magnitude at 
this point [11111110. The large aspect ratio of CNTs 
means that, for a disordered distribution of positions and 
orientations, a single particle provides connectivity over a 
longer distance than several more spherical particles with 
the same total volume. Hence, an insulating matrix can 
be loaded with a smaller volume fraction of conducting 
filler, leading to a lighter material. 

The percolation threshold of the CNTs depends both 
on their aspect ratio, and on the interactions between 
them. Recent studies show that the addition of sur- 
factant to a suspension of CNT bundles (CNTBs) sub- 
stantially lowers the percolation threshold 6]. At suffi- 
ciently high concentrations, the surfactant self-assembles 
into micelles, which have been shown to induce depletion 
attraction between the bundles Q . Due to the weakness 
of the attraction, the percolation is reversible and the 
system remains at thermodynamic equilibrium, in con- 
trast to strongly interacting particles, which may form 
structurally arrested gels [8|. 

The depletion forces arise from competition between 
the translational and orientational entropy of the CNTBs 
on the one hand, and the translational entropy of the mi- 
celles on the other. If two bundles approach each other 
closely, the overall volume available to the micelles in- 
creases, thereby inducing an effective attraction. For 
rod-shaped particles like CNTBs, this depletion effect is 
highly anisotropic, since the total volume excluded to 
the depletant depends on the relative orientations of the 
bundles 0, ^| . The excluded volume is minimized by 



parallel arrangements of neighboring rods, but alignment 
also reduces the spatial extent of a cluster, and should 
therefore suppress rather than enhance the formation of 
percolating networks [TH ]. Even in systems of spheres, it 
is possible for attraction to raise the percolation thresh- 
old Nevertheless, recent analytic work on mixtures 
of rods and spherical depletants predicts that depletion 
does lower the percolation threshold of the rods for all 
depletant concentrations [l^j ]. In the present study, we 
will see that depletion induces relatively local changes in 
the structure of an isotropic suspension of rods, but that 
these changes significantly affect the extent of clustering 
and therefore the position of the percolation threshold. 

We model the CNTBs as hard spherocylinders, i.e., 
cylinders of length L and diameter D capped by hemi- 
spheres. To make as few assumptions as possible about 
the effect of depletion, we choose to model the surfactant 
micelles explicitly rather than by employing an effective 
potential. Explicit representation has the advantage of 
automatically including many-body effects, due to mul- 
tiple exclusion of the same volume, which can never be 
perfectly accounted for by a pair potential. In the exper- 
iments, the diameters of CNTBs and micelles are com- 
parable 0, 0, and we therefore model the micelles as 
hard spheres of diameter D. The range of the depletion 
interaction does depend on the depletant size, but for a 
given sphere volume fraction, a spot check at aspect ratio 
L/D =10 revealed only a 2% change in the percolation 
threshold using sphere diameters of 0.5D to 1.3D. 

Similar mixtures have been the subject of other theo- 
retical and simulation studies [jj, 0, [y| . In particular, 
the phase diagram for low volume fractions and roughly 
equal rod and sphere diameters was studied within the 
second virial approximation by Matsuyama and Kato 
[lfij . Most other studies, however, have been concerned 
with different regimes, either where the rods are in a liq- 
uid crystalline, rather than isotropic phase |l6| , or where 
the rods are re gard ed as the depletant rather than the 
spheres ^3, ITsl Il9l| . We believe the present work to be 
the first time that depletion-induced percolation has been 
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simulated in rod-sphere mixtures. 

Monte Carlo simulations were performed at fixed vol- 
ume and number of particles using displacements and 
(for rods) rotations of individual particles. The hard po- 
tential reduces the Monte Carlo acceptance criterion to 
the rejection of any trial move that generates an overlap 
between rods and/or spheres. For efficient overlap detec- 
tion, the simulation box was divided into cells with edge 
length I > D |2jj. A cubic simulation box with periodic 
boundary conditions was used and, except where stated, 
the box length was 5L. A typical simulation involved 
5 x 10 5 trial moves per particle (sweeps) for equilibration 
and a further 5 x 10 5 for sampling. Overall alignment 
of the rods was monitored using the nematic order pa- 
rameter S, which is the largest eigenvalue of the tensor 
Q = (2N I )~ 1 J2f' (3uiUi — I). Here, is a unit vec- 
tor along the axis of rod i, I is the identity matrix, and 
the sum is over all N T rods in the system. S is unity 
in a fully aligned configuration and zero in a perfectly 
isotropic fluid. 

To define clusters of rods, a connectivity criterion is 
required. We are not concerned with details of electrical 
conductivity at CNTB contacts, and regard two rods to 
be connected if their surfaces approach closer than 0.2D, 
which is equivalent to the line segments at the sphcro- 
cylinders' axes approaching closer than A = 1.2D. This 
is a physically reasonable, though necessarily arbitrary 
cutoff. A cluster was counted as percolating if its parti- 
cles were connected to their own periodic images. Perco- 
lation of the configuration was tested every 20 sweeps. 

The aspect ratio of CNTBs can exceed 100, which is 
beyond the reach of our simulations. However, the emer- 
gence of long rod behaviour in the rod-sphere mixture be- 
comes clear by simulating increasing L/D starting from 
relatively small values. Figure demonstrates the ef- 
fect of depletion forces on the percolation of rods with 
L/D = 6. The equilibrium probability P(t) t ) of finding 
a percolating cluster is plotted as a function of the rod 
volume fraction rj r = N T Vh,/V, where N T is the number 
of rods, i> h = 7r£> 3 (2 + 3L/D)/12 is the volume of the 
rod's hard core, and V is the volume of the simulation 
cell. Addition of a volume fraction i] s = 0.1 of spheres 
shifts P(rj T ) to lower rod volume fraction by 17.5%. 

In an infinitely large system, P{rj T ) would rise instanta- 
neously at the percolation threshold, but the simulations 
measure a sigmoidal curve whose width decreases with 
increasing cell size due to finite size effects, as demon- 
strated for the pure rod case in Fig.^ Nevertheless, the 
volume fraction at which P(rj T ) passes through 0.5 is al- 
most independent of the cell size and we therefore adopt 
P(v? cr ) — 0-5 as a robust definition of the percolation 
threshold rj^ ei . In general, above the sharp onset of per- 
colation at 7yf er , the conductivity of a network continues 
to increase with critical exponent t jUI^ as (r? r — »?? er )*- 

As expected, r]f eT is lower for longer rods, as shown 
for a selection of L/D ratios in Table [I] Data are given 




volume fraction of rods, T\ 

FIG. 1: Solid lines: percolation probability as a function of 
rod volume fraction for L/D = 6 and two sphere volume frac- 
tions rj s in a box of edge 5L — 30D. Dashed and dotted lines: 
equivalent results for rj s = only at box lengths 2073 and 
45D, respectively. A typical uncertainty in the probability 
is 0.004. Inset: percolation threshold as a function of sphere 
volume fraction for L/D = 6. 



TABLE I: Percolation threshold for selected aspect ratios 
L/D at two sphere volume fractions r] s . rif BT is the volume 
fraction of rods at percolation, and n pcr is the average num- 
ber of (rod) neighbors per rod. S is the fractional decrease in 
r/f eT caused by sphere volume fraction ?7 S = 0.03. 
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for rods with no depletant, and with a volume fraction 
77s = 0.03 of spheres, which is comparable to the highest 
concentration of micelles used in Ref. 0- The fractional 
decrease S in the percolation threshold due to the spheres 
is 5.2% at L/D = 4, rising gradually to 8.3% by L/D = 
25 in contrast to theoretical predictions of constant 5 [l3j |. 
The increasing shift is likely to continue as the aspect 
ratio rises towards that typical of CNTBs. However, if a 
scaling regime of the form r]f cr oc (L/D)~ a is entered for 
sufficiently long rods, and the value of a is not affected 
by the presence of spheres, then the fractional drop in 
7]f eI caused by a given volume fraction of spheres would 
become independent of the aspect ratio. 

In fact, such a scaling has been predicted for ideal sys- 
tems of freely interpenetrable rods. In the hope of esti- 
mating percolation thresholds without the need for de- 
tailed calculation, several quantities have been suggested 
as approximate dimensional invariants at the threshold. 
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Candidates included the number of bonds per particle 
and the total excluded volume [3, 0] , and it has been 
shown that the zero of the Euler characteristic correlates 
well with percolation [2H l2t| . Although none of these 
approaches works universally well for arbitrary particle 
shapes and mixtures, the quantity p pcr (v cx ) is an approx- 
imate invariant for ideal rods [2^] , where p per is the num- 
ber density at percolation and (v ex ) is the so-called ex- 
cluded volume (orientationally averaged) that would be 
inaccessible to the center of one particle due to the pres- 
ence of another if overlaps were forbidden. Defining the 
total critical volume by V c = p pcl v, where v is the volume 
of an ideal rod, the invariant becomes V c (f ex ) /v. For long 
rods, (v cx )/v oc L and hence V c oc l/L. 

Ideal rods rapidly enter this scaling regime [22^. as 
shown by the dotted line in Fig. [2J Each evaluation of 
the percolation probability was obtained from 10 000 un- 
corrected and randomly generated configurations of line 
segments, using the same box size and connectivity cri- 
terion as in the spherocylinder simulations. The intro- 
duction of a hard core induces local structure into the 
isotropic fluid and rules out a large fraction of percolat- 
ing ideal configurations, raising the percolation thresh- 
old. The total critical volume for hard rods with and 
without depletant spheres is also shown in Fig.|3 For the 
aspect ratios accessible to our simulations, no power-law 
scaling is observed. However, theoretical work [27| pre- 
dicts that the effect of a hard core is to postpone scaling 
to larger aspect ratios and that 1 /L scaling is eventually 
recovered. Table [i] shows a decline in the rate of increase 
of S with L/D, and in the scaling regime S would level 
out to a constant. A linear extrapolation of 5 to experi- 
mental aspect ratios starting at L/D = 200 therefore 
provides a crude upper bound to the contribution of pure 
entropy in the CNTB-micellc experiments. The extrapo- 
lation predicts a lowering of 30-40%, a sizeable effect but 
far smaller than the experimental one. The additional 
influence of direct energetic interactions in the CNTB 
suspension could be tested by replacing the micelles by a 
different depletion agent with the same excluded volume. 

We now turn to structural characteristics of the fluid. 
In all the simulations the rods remain in the isotropic 
phase, which has no long-range order in either the center- 
of-mass positions or the alignment of the rods. However, 
the hard interactions give rise to local structure, which 
is affected by the addition of hard spheres. Local align- 
ment is revealed by the orientational correlation function 
92{r) = (P 2 (u* • Uf))r, where P 2 {x) = (3x 2 - l)/2 is 
the second Legendre polynomial, and the averaging is re- 
stricted to pairs of rods whose centers are separated by 
Tij = r. As shown in Fig.|3Ja), gi (r) decays to zero in an 
isotropic fluid, but becomes positive when rods approach 
closely enough to restrict each other's rotation. 

The addition of hard spheres to the rods encourages 
nearby rods to align, since alignment decreases the to- 
tal volume excluded to the spheres by an amount pro- 
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cylinder aspect ratio L/D 

FIG. 2: Dependence of the total critical volume on the aspect 
ratio of the rods for sphere volume fractions 77,5 = (solid line) 
and rj s — 0.03 (dashed line), and for the ideal system of freely 
interpenetrable rods (dotted line). 



portional to the length of the adjacent segment, thereby 
increasing the entropy of the spheres at the expense of 
some orientational entropy of the rods. Figure|3fa) shows 
that the spheres greatly enhance alignment of rods with 
nearby centers and (inset) that the centers are pushed 
closer together. The range of the effect is on the order of 
the sphere diameter D and decays to the underlying cor- 
relation, whose range depends on L. For sphere concen- 
trations = 0.03 around those of the micelles in Ref. 6, 
there is very little increase in alignment beyond the range 
D, and no significant bundling of sphcrocylinders was 
observed. At higher rj r , the spread of depletion- induced 
alignment lowers the bulk isotropic-nematic transition 
density 0], but following the percolation line r] pcl '(r] s ) 
(inset of Fig. , we find that the nematic is preempted 
by rod-sphere demixing (around ?7 S = 0.22 for L/D = 6). 

To understand cluster formation and percolation, a 
crucial structural property is the probability distribution 
p(d) of shortest distances d between pairs of rod axes, 
since it is the criterion d < A that defines contact be- 
tween two rods. To reveal the effect of the hard poten- 
tial on the structure, we divide p{d) by the equivalent 
distribution Pid{d) in a system of ideal rods. The ra- 
tio p(d)/pia(d) is plotted in Fig. IHIb). As with the first 
peak in the radial distribution function of a simple fluid, 
the distribution of axis-axis separations rises above the 
random distribution for small d. The addition of hard 
spheres again has a large but short-ranged effect on the 
distribution, pushing the rods together at any point of 
approach along their lengths. For the high sphere vol- 
ume fraction of ?7 S = 0.15 a small second peak emerges 
at d « 2. 2D, indicating a weak correlation between two 
contacts with an intervening rod or sphere. The enhance- 
ment of p(d) by the spheres for separations d < A is a 
highly localized adjustment to the structure of the pure 
rod fluid, but increases the connectivity enough to cause 
the cluster size to diverge at lower rod concentration. 
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FIG. 3: Local structure of a spherocy Under fluid with 
L/D = 6 at volume fraction r] T = 0.142. (a) Orientational 
correlation function g2(r), inset: radial distribution function, 
g(r), (b) distribution p(d) of axis-axis separations normal- 
ized by the ideal distribution Pi<x(d). Results for three sphere 
volume fractions are shown: r/ s — (solid lines), 7/ s = 0.1 
(dashed) and 7/ s = 0.15 (dotted). 

The average number n pcr of contacts per rod at per- 
colation (Table P), decreases with the aspect ratio, con- 
firming that, as for ideal objects, this quantity is not an 
invariant [221 1241 |25j . The number of contacts is a sen- 
sitive function of the rod density, so it is interesting to 
note that n por has similar values at 77 s = and 0.03; the 
spheres lower the percolation threshold, but they increase 
the number of contacts per rod at the lowered threshold 
roughly to its percolation value their absence. The low 
n per for long rods show how sparsely connected the net- 
work is at the onset of percolation. The strength of the 
network continues to increase sharply above ?y per . 

Using a simple model of hard spherocylinders and 
spheres, we have isolated the effect of depletion on the 
percolation threshold of nanorod networks. The spheres 
induce pronounced but short-ranged changes in the struc- 
ture of the isotropic rod fluid, including enhanced local 
alignment. At depletant concentrations relevant to ex- 
perimental work on suspensions of carbon nanotubes and 
micelles Q , the alignment is too weak to cause global ne- 
matic order or significant aggregation of rods. Although 
pairs of particles extend less far when mutually aligned, 



depletion enhances contacts at all positions along the 
rods, leading to larger clusters and a decrease in the per- 
colation threshold. It should therefore be possible to en- 
hance the percolation of nanorods using depletants other 
than micelles, such as polymer coils or latex spheres. The 
fractional decrease caused by a given volume fraction of 
depletant gradually increases with the aspect ratio of the 
rods. This prediction could, in principle, be tested if 
shorter nanorods can be synthesized and characterized. 
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